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We study the dynamics of a micromaser where the pumping atoms are strongly driven by a 
resonant classical field during their transit through the cavity mode. We derive a master equation 
for this strongly-driven micromaser, involving the contributions of the unitary atom-field interactions 
and the dissipative effects of a thermal bath. We find analytical solutions for the temporal evolution 
and the steady-state of this system by means of phase-space techniques, providing an unusual 
solvable model of an open quantum system, including pumping and decoherence. We derive closed 
expressions for all relevant expectation values, describing the statistics of the cavity field and the 
detected atomic levels. The transient regime shows the build-up of mixtures of mesoscopic fields 
evolving towards a superpoissonian steady-state field that, nevertheless, yields atomic correlations 
that exhibit stronger nonclassical features than the conventional micromaser. 

PACS numbers: 42.50.Pq, 03.65.Yz, 32.80.Qk 



I. INTRODUCTION 

The micromaser P| is a fundamental system in cav- 
ity quantum electrodynamics (CQED) (J| where a sin- 
gle cavity mode is pumped by a poissonian beam of ex- 
cited two-level atoms, in such a way that at most one 
atom is inside the cavity. While each atom traverses the 
cavity, the entangled atom-cavity system undergoes the 
Jaynes-Cummings (JC) interaction Q and, when this 
is not the case, the field decays due to its coupling to 
the environment. It has been demonstrated that the mi- 
cromaser exhibits interesting nonclassical features in the 
field statistics and in the atomic correlations, like sub- 
poissonian photon statistics trapping states 's^ and 
number states i^] of the cavity field, and antibunched 
atomic correlations 0- 

Recently, it was shown that the effective coupling be- 
tween an atom and a single cavity mode can be drasti- 
cally modified in the presence of a strong external driving 
field 'il. For example, under resonant conditions among 
the atomic transition, the cavity mode and the external 
field, it is possible to engineer a resonant JC and anti- 
JC interaction simultaneously. In this case, the usual 
atom-field Rabi oscillations do not play any further role, 
leaving their place to conditional field displacements de- 
pending on the atomic internal states: "Schrodinger cat 
states." Generation of entanglement and its measure- 
ment is a very active field in CQED with diverse implica- 
tions in fundamental and applied quantum physics 0, Q . 
From that point of view, a more realistic study of this 
new scheme and its different variants, where dissipative 
effects are fully taken into account, is desirable. 

In this paper we develop the theory of a strongly-driven 
micromaser (SDM) and study the new features of the 
field statistics and the atomic correlations. We formulate 
a master equation for the cavity field density operator 
where the gain originates in the coherent interaction of 



the field with two-level atoms which are strongly driven 
by an external driving while they are inside the cavity . 
The losses in the master equation stem from the interac- 
tion of the cavity field with a thermal bath. The addi- 
tional strong driving acting on the atoms changes dras- 
tically the SDM dynamics as compared with the conven- 
tional micromaser. 

We show that the SDM master equation can be solved 
analytically by means of phase-space methods, providing 
an unusual solvable model of an open quantum system 
under realistic conditions. In this way we are able to 
trace the temporal evolution of the SDM field from an 
arbitrary initial state to the final steady-state, which hap- 
pens to be superpoissonian. We derive closed expressions 
for the main quantities which characterize the statistics 
of the cavity photons and the detected atoms. We find 
that, despite the classicality of the SDM field steady- 
state, the atomic correlations exhibit stronger nonclassi- 
cal features when compared with the conventional micro- 
maser. The description of system dynamics is illustrated 
also by numerical results which support and complement 
the analytical ones. It is worth noting that the present 
SDM system differs markedly from the recently investi- 
gated coherently-driven micromaser where an ex- 
ternal field continuously drives the cavity mode and no 
strong driving regime is required. 

In Sec. we present the atom-field Hamiltonian of 
the SDM and illustrate its main dynamical features. In 
particular, we show how the cavity field changes as a 
consequence of the passage of one, two, or more atoms. 
We pay careful attention to the crucial difference between 
the situations in which the atoms remain unobserved and 
in which their final state is detected. Then, in Sec. lIIII we 
formulate the SDM master equation and derive its time- 
dependent analytical solution. As an application, we use 
it to calculate the expectation values of various important 
observables of the photon field, and of the atom counting 
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statistics. Numerical results are discussed in Sec . 11 VI and 
conclusions are drawn in Sec. ]V\ 

II. HAMILTONIAN AND COHERENT 
DYNAMICS 

The fully resonant interaction between one mode of a 
high-Q cavity and a two-level atom strongly driven by a 
classical external field can be described by the following 
Hamiltonian |3| 

H^^ia^+a)ia^+a), (1) 

where g is the atom-cavity mode coupling constant, a (a^) 
the field annihilation (creation) operator, and a = |5)(e| 
(ct^ = |e)((ji|) the atomic lowering (raising) operator. The 
Hamiltonian of Eq. written in the interaction picture, 
was derived in 8] in the strong driving regime f2 3> g, 
where D, is the Rabi frequency associated with the exter- 
nal field. 

A first feature of Hamiltonian in Eq. is the presence 
of resonant terms of both the Jaynes-Cummings type, 
a^a + cra'l', and of the anti-JC type, a^a^ + aa. Actu- 
ally, the latter terms are usually negligible in CQED, 
whereas they can be of importance in other systems, like 
ion traps A second feature, that is discussed below, 
is the direct generation of the so-called Schrodinger cat 
states of the cavity field. 

The Hamiltonian of Eq. ^ generates a time evolution 
that is described by the unitary operator 

u{o^m\+){+\+Di-o\-){-\, (2) 

where the parameter ^ = —igT/2 is proportional to the 
interaction time r; D{a) = exp {aa^ — a*a) is the uni- 
tary displacement operator; and |±) — {l/y/2){\g) ± |e)) 
are the eigenstates of the atomic operator ax = cr^ + cr 
with eigenvalues ±1, respectively. Starting from an ini- 
tial state with the cavity field in a state pF,o and one 
atom injected in the excited state, 

Po = PF,o «) |e)(e| , (3) 

the density operator evolves to the state 

pi = U{OpoU{-0 ■ (4) 

We consider two cases, one in which the atom is not 
observed when it leaves the cavity and the other in which 
it is detected in the upper or lower level. 

A. Atoms are not observed 

If we are interested in the cavity field after the atomic 
transit without measuring the state of the outgoing atom, 
the cavity field density operator is 

PF,i = TvAPi = i [D{Opf,oD{-0 + D{-OpF,oDiO] , 

(5) 




FIG. 1: Wigner distributions Wi{a) and W2{a) of the cavity 
field, see Eq. @, in the case of one and two crossing atoms, 
respectively, provided that no attempt is made to determine 
the final internal state of the atoms. Here, as well as in all 
the following figures, |^| = tt, the cavity field is initially in the 
vacuum state, and the atoms in the excited state. 

where Tr^i is the partial trace over the atomic degrees 
of freedom, and we used Eqs. Q-Q. If the cavity is 
initially in the vacuum state, pF,o = |0)(0|j its state turns 
into 

PF,i = i(io(ei + i-o(-ei) , (6) 

which is a statistical mixture of two coherent states with 
the same amplitude and opposite phases. 

It is of interest for what follows to describe the time 
evolution of the cavity field in phase space. While for 
the vacuum state the Wigner function is the Gaussian 
Wo{a) = 2exp(— 2|ap), for the state of Eq. ^ we have 

Wi{a) = 2exp [-2(|ap + I^H] cosh(4|C|Im a) . (7) 

This last expression shows that the interaction with a 
driven atom can affect the rotational symmetry of the 
initial vacuum state. The section of Wi(q!) along the 
real axis remains a Gaussian with the initial minimum 
uncertainty, while the section along the imaginary axis 
is broadened, showing a two-peaked structure that re- 
sults from the superposition of two Gaussians centered 
at a = The two peaks are well resolved if |^| » 1/2, 
which requires that the system is operated in the strong- 
coupling regime (gr ^1). 

If the cavity mode, initially in the vacuum state, inter- 
acts with 1, 2, . . . , m strongly driven atoms, all assumed 
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to pass through the cavity (one by one) with the same 
interaction time r and to leave it unobserved, the final 
field state is given by the density operator 



— Y 

n=0 



|(m- 2n)0((m- 2n)^| , (8) 



which describes a mixture of ?7i + 1 coherent states. This 
state is represented in phase space by the Wigner func- 
tion 



Wmia) 



1 



2m — 1 / ^ 

71 = 



^ exp [-2 



\a-{m- 2n)^|^ 



(9) 

Hence, in the strong coupling regime and with neg- 
ligible dissipation effects, multi-peaked cavity field dis- 
tributions can be generated in phase-space, where the 
number of peaks increases with the number of driven 
atoms injected in the cavity. The subsequent genera- 
tion of such cavity field states is illustrated in Fig. ^ 
where we show the Wigner functions Wi (a) , W2 (a) cal- 
culated from Eq. 0. After the transit of m atoms the 
function Wm{<x) exhibits m + 1 peaks, all centered on 
the imaginary axis and with a center-to-center distance 
of 2|^| — gr. If m is even, the peaks are centered at 
Rea — and Ima = 0, ±2|^|, . . . , where the central peak 
in the origin is the highest one. If m is odd, the centers 
are at Rca ~ and Ima = ±|^|,±3|f|, . . ., the highest 
peaks being at a = 



B. Atoms are detected 

Now we consider the cavity field in the case in which 
the strongly driven atoms are detected when they leave 
the cavity, e.g. their state is determined by selective field 
ionization, where for simplicity we assume perfect detec- 
tor efficiency. We are again interested in the cavity field 
properties after interaction. Proceeding again from the 
initial condition of Eq. and using Eqs. (Hjl-Q, the 
cavity field density operators after detecting the atom in 
the excited state p'p\ or in the ground state pf,\ are 



(e) 

Pfi 



ig) 



2 [Rex(20 + 1] 

[D{0-D{-e}]PF.o [D{0-D{~0] 
2 [Rcx(2C) - 1] 



(10) 



where x{P) = Tr_F {pFfiD{fi)} is the characteristic func- 
tion for symmetrical ordering of the field operators 
and Tvp is the partial trace over the field variables. If 
the cavity is initially in the vacuum state, then 

(e) _ |0(el + |-0(-el + IO(-el + |-0(CI 

Pfs — 



(9) _ IO(ei + 

Pf,i — 



2[l + exp(-2|^|2)] 
2[l-exp(-2|eP)] 



which are pure states of the cavity field instead of the 
statistical mixtures of Eq. © . Actually, the cavity field 

state vectors are IV'-F.i)*'^'''''^^ oc (|^) ± (| — C)); that is 
the superposition of two coherent states of the kind gen- 
erated and monitored in the dispersive regime of cavity 
QED in . More elaborated superposition states were 
investigated in Ref. Q , where also many-atom states were 
considered. The mean photon number of the field states 
of Eq. lfTT|l are 



(^1 



lTexp(^2|CP) 
l±exp(-2|C|2) 



(12) 



whereas (A^i) = in the case of an unmeasured atom 

of Sec. urn 

The Wigner functions representing the states of 
Eq. lfTT|l are 



cosh(4|^|ImQ;) 



±cos(4|^|Rea)]/(l±e 



-2|?|- 



(13) 



Beyond the two-peaked structure, present in Eq. |(7J) for 
an unmeasured atom (see also Fig. Pl, the presence of 
the sinusoidal interference term implies that the Wigner 
functions in Eq. H13|l can exhibit strong oscillations with 
period 7r/2|^|. They can even take negative values (see 
Fig. [2J), which is a signature of the quantum nature of 
the cavity field states of Eq. (|ll|l . In particular, in the 

origin of phase space, Wi'^^'^^\o) — ±2. 

If a second atom crosses the cavity and is also detected 
in the upper or lower state soon after, the cavity field is 
projected onto one of the following states. 



(ee) 
Pf,2 



Pf,2 



Agg) 

Pf,2 



[\20 (2^1 + 4|0) (0| + I - 20 (-2^1 + \20 (-2^1 
+ |-2e)(2e|+2(|2O(0| + |0)(2e|-l-|0)(-2C| 
+ |-2e)(0|)]/[2(3 + 4e-2|«lVe-8|«l^)], 

= [120(2^1 + 1-20 (-2^1 
-|20(-2e|-|-2e)(2e|]/[2(l-e-«l«l^)], 

[120 (20 + 4|0) (0| + I - 20 (-20 + |20 (-20 
+ 1 - 20(20 - 2(|2O(0| + |0)(20 + |0)(-20 



+ |-20(0|)]/[2(3-4e 



-2|?|- 



(14) 



(11) 



where, for instance, (eg) means that the first atom is de- 
tected in the upper state and the second one in the lower 
state. The expressions of Eq. (|14|) should be compared 
with the statistical mixture of Eq. ((SJ in the case of two 
unobserved atoms (m = 2). The density operators in 
Eq. (|14(l describe mesoscopic Schrodinger-cat-like states 
of the cavity field, the corresponding state vectors being 

|^^^^)(ee),to) ^ (|2O±2|0) + |-2O), 

|V'W^'^^ = |V'K2)'^'^(x (120-1 -20). (15) 
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(a) 



(a) 




■8 -4 "'^ Rea 



(b) 




FIG. 2: Wigner distribution of the cavity field for the case of 
one crossing atom detected (a) in the excited state, W^''^'(a); 
(b) in the ground state, W^^\a); see Eq. 11311 . 



The Wigner functions which represent the above states 
in phase space, whose behavior is depicted in Fig. |3| can 

be written as 



W^^'^a) = W^^'^a) = 2e-2|"l'[~cos(8|e|Rea) 



2e-2|"l' [2 + e-**l«l' cosh(8|^|Im a) 
+4e"2|«l'cosh(4|^|Ima)cos(4|^|Rea) 
+ cos(8|C|Re a)] / (3 + 4e~2|?|^ _^ p-8|«|^ 



rise 



-e"^l«l cost 



3h(8|e|Inia)]/(l 
2e-2|"l' [2 + e"^l«l' cosh(8|^|Ima) 
-4e"2|«l' cosh(4|^|Ima) cos(4|^|Rea) 



3(8|C|Rea)]/(3-4e-2|?l^ 



(16) 

We consider now the atomic statistics at the exit of 
the cavity independent of the cavity field state. If an 
atom is injected in the cavity, so that the initial state is 
Pf ^ |e)(e|, the atomic state after the interaction will be 

PA = TTF{U{OpF<E>\e){e\U{-0} 
= i[u + Rex(20(|e)(e|-|5)(.g|) 



+zlmx(2e)(|<?)(e|-|e)(<?|) 



(17) 




(b) 



(c) 



Ima 





I met 



-10 -2 



Re™ 




-10 -2 "Rea 




-10 



Rea 



FIG. 3: Wigner distribution of the cavity field for the 
case of two crossing atoms detected (a) both in the excited 
state, W^'"'\a); (b) one in the excited state and the other 
in the ground state, W2'^^\a); (c) both in the ground state, 
W^^^\a)- see Eq. llT^ . 



Here, Tip is the partial trace over the cavity mode de- 
grees of freedom. In the derivation of Eq. H17(l we used 
also Eqs. Q and Q as well as the property x(/3) = 
x(— /?)*. Finally, we can calculate the probability Pe,g 
for atomic detection in the upper or lower state as 



Pe,g = {e,g\pA\e,g) = - [l ± Rex(2e)] • (18) 
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III. MASTER EQUATION AND ANALYTICAL 
RESULTS 



A. SDM master equation and analytical solutions 

We consider now the case of a poissonian beam of two- 
level Rydberg atoms, with pumping rate r, interacting 
with a single mode of a microwave high-Q cavity. While 
the atoms are inside the cavity, during an interaction time 
T, they are strongly driven by an additional classical field, 
the whole system being in resonance. We study a regime 
where at most one atom is present inside the cavity, that 
is, r r~^, and where the decay of atomic Rydberg 
levels is negligible. Between two successive atoms, the 
cavity field decays due to its interaction with a thermal 
bath. Hence, even though the present strongly-driven mi- 
cromaser (SDM) looks quite similar to the conventional 
micromaser, there is a crucial difference in the pump- 
ing dynamics. As we discussed in the previous section, 
the JC interaction and its Rabi oscillations do not rule 
any more the unitary atom-field interaction. Their role 
is taken over by field displacements conditioned on the 
atomic internal states. 

The dynamics of the SDM field is ruled by the interplay 
between the amplification process, due to the interaction 
with the driven atoms, and the dissipation process, occur- 
ring when no atom crosses the cavity. In a coarse-grained 
description the gain rate of the cavity mode follows from 
Eq.®, 



dpi 



dt 



= - [D{^)pfD{-0 + D{-£,)pfD{^) - 2pf\ 



gain 



(19) 

The loss rate, due to the interaction of the cavity field 
with a thermal bath , is given by 



dph 



dt 



loss 



7(» + 1) r t o t , t 1 

^ [a'apF — lappa + ppa cl\ 

--^ \aa) pf — 2a^ ppa + pi?aa^] = Cpi. 



(20) 

where 7 is the cavity photon decay rate and n the mean 
thermal photon number. By combining Eqs. Ijiyfl and 
(|20|l we obtain the SDM master equation 



dpp 
dt 



-[D{OPFD[-i) + D{-OppD{0 - 2pp] + Cpp . 

(21) 

At variance with the conventional micromaser, whose 
master equation has only an analytical solution in steady- 
state the SDM master equation (|2f|l can be fully 
solved for any time. The key step is a transformation 
of the equation of motion H2fl) for the density operator 
into the corresponding partial differential equation for its 
symmetrically ordered characteristic function x(/3), al- 
ready introduced in the previous section. We recall that 



x(/3) is related to the field density operator by 

pp^n-' J XmD~\f3)d'p (22) 

and that the Wigner function, W{a), is the 2D Fourier 
transform of the characteristic function y(0) 12\. 

Using known operator techniques |fj|, we thus map 
Eq. H2f |l into the following partial differential equation 
for x{f3,f3*,t), 



d_ 

dt 



(2n-M)|/3|2 



x(A/3*,i) 



| + /3*|,)x(/?,r,t). 



(23) 



We transform Eq. (|23|l to cartesian coordinates by setting 
^*/3 = 1^1(2: -f iy), which gives 



d 7 
dt^2 



where 



d d i\ 
'^-Q^ + y-Q^ + a;G"(a;) + yF'{y)\ j x(x, y, i) = , 

^ (24) 



G{x) 

ny) 



2nb + 1 o 

X 



AN,, Sin2(|e|z) + (2nfc + f)z2 



dz , (25) 



and Nex = r/7. The solution of Eq. 1)24(1 reads 



x{x,y,t) 



x''''''Hx,y)xo{xe ^2 ^ye I) 



Here, 



(26) 



(27) 



is the steady-state solution of Eq. H24|) that is reached 
for t 00, and xo = T^t^p[pp{0)D] is the characteristic 
function corresponding to the initial field state pp{0)- By 
substituting the expressions of Eq. H25() into Eq. (|27|l . we 
arrive at the steady-state characteristic function 

X^"'Hx,y) = exp(-(n + ^){x^ + y^) 

+4iVe4-7e-ln(2|^|2/) + Ci(2|C|y)]) , 



(28) 

where 7e is Euler's constant and Ci(2|^|?/) is the Cosine 
Integral defined as 



my 



Ci(2|e|?/) = 7e + ln(2|e|y) 



cos(z) — f 



dz . (29) 



6 



Note that starting from a real initial function, the time- 
dependent solution in Eq. H2()|) will always produce a 
real characteristic function. This follows from the in- 
variance of the time evolution equation in Eq. (|23|l un- 
der the transformation P —^ — /3, and from the property 
xiP) = x(-/3)*- Accordingly, the solution of Eq. if^ 
depends only on the modulus of the (imaginary) param- 
eter ^, which is in agreement with the invariance of the 
master equation (|21|l under the transformation ^ — ^. 



B. SDM Held statistics 

It is clear that from the general solution of the charac- 
teristic function, Eq. (|26|l . we can calculate the field den- 
sity operator prit), following Eq. (|^ . and have access 
to the field statistics. However, this information can be 
directly extracted from the characteristic function itself, 
since the expectation value of any symmetrically ordered 
product of the operators a and is given by 12] 



ryn^n 



9'"+"x(/3,/3*) 



'sym 



(30) 



f3=l3'=0 



For example, we have 



(ait)) 
{a\t)a{t)) 
(a^t)) 



(a(0))e-^*, 

(at(0)a(0))e-^* + (7Ve,|CP 



n)(l 



(a2(0))e-^*+iV,,e'(l-e-^*). 



(31) 



From the expressions in Eq. H31() and their complex 
conjugates, we derive the steady-state expectation val- 
ues of the field amplitude, photon number, and quadra- 
ture variances (Aa;i_2)^, with xi — (l/2)(a^ + a) and 
X2 = {i/2){a^ ~a), 



V (ss) 



[Axi^'^'^f = i(l 

{AX2 



(ss)n2 



0, 

2n), 

i(l + 2n + 4iVex|^n 



(32) 



We see that at steady-state the expectation value of the 
SDM field is zero, whereas the mean photon number is a 
quadratic function of the modulus of the single atom dis- 
placement parameter ^. Also, for vanishing cavity tem- 
perature (n — > 0) , the variance of the quadrature oper- 
ator xi remains at the minimum value 1/4, whereas the 
variance of the orthogonal quadrature X2, the one that 
is being driven by the system, is broadened by a factor 
equal to (a^a)^^'^\ 

Another important quantity describing the photon 
statistics is the Fano-Mandel parameter Q = [(AA^)'^ — 
{N)]/{N), where N = {a'' a). For the SDM at steady- 
state we obtain 



Q = l + iVe,|e|'+: 



(33) 



which describes a persistent superpoissonian behavior 
(Q > 1). This result is at variance with the conventional 
micromaser steady-state which alternates between 
superpoissonian and subpoissonian statistics. 



C. Atomic correlations 

Just as it is the case for the micromaser, in the SDM we 
are not able to measure the cavity field directly. There- 
fore, we use the atoms not only for pumping the cav- 
ity mode but also as a source of information about the 
SDM field. The theory of the detector clicks statistics as 
well as the connection between the statistics of the de- 
tected atoms and the cavity field was developed for the 
conventional micromaser in Ref. jl^ . According to this 
approach, the detection of the exiting atom, initially ex- 
cited, in the ground state is described by an operator A 
whose action on the field density operator pF is 

ApF = ^[Di-OPFDiO + DiOpFD{~0 

-Di-OpFDi-0 - DiOpFDiO] . (34) 

Likewise, the click operator B for detection in the excited 
state acts as follows 

1 



BpF = - [Di-OPFD{0 + DiOpFDi-0 

+d{-OpfD{-0 + d{OpfD{0] ■ 



(35) 



Note that the above expressions have the same structure 
as the field state operators of Eqs. (|10|l . The probability 
to detect the atom in the ground (excited) state Pg (pe) 
can be calculated as Ttf[Apf] {T^^FiBpF]), giving the 
results of Eq. H18|) , which can be further simplified after 
the derivation of a real expression for the characteristic 
function x, so that 



(36) 



The introduction of the above click operators enables us 
to describe the atomic correlation functions, or condi- 
tional probabilities for two consecutive detector clicks. 
For instance, the correlation function for a detection of 
the first atom in the excited state after a detection of 
the second atom in the ground state separated by a time 
interval t is given by 



1 + X{2U) 



TrF[6p("")]Tri.[yip(^^)] 1 + X^'^'^H'^O 



(37) 



where Cq — C + r{A + B — 1), x^^^^ is the steady-state 
characteristic function |2Hll, and 



2x(""-'(2Ce"^*) - 1 - x^'"-'(4Ce"^*) 
2_2x("'^)(2^e-4t) 

xexp(-(n+i)|2Cp(l-e-4*)). (38) 
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FIG. 4: Atomic correlations Geg with parameters N^^ ~ 50, 
n = 0.03. The figures are plotted for time intervals t — 
(top) and t — 0.01 (bottom) between detection clicks. 



The two-click correlation function in Eq. H37|l is the ra- 
tio between the conditional probability to have a second 
e-click at time t after a first g-click occurred at time 0, 
and the probability for an e-click when the cavity field 
is in the steady-state. All other correlation functions for 
click pairs Ge 
expressions 



Ggg and Geg are given by the following 



Gee{t) 



Trf [^e^t'Up^'^'^)] _ J_ - x(2?, t) 
Tri.[y^p(^'')]Trj.[Sp(^^)] 



where 



(39) 



2 + 2x(""H2?e-^*) 



X exp 



(-(n+i)|2^|2(l-e'4*)). (40) 



However, like in the case of the conventional micro- 
maser, the two-click correlation functions of Eqs. 1)3 7|) 
and (|39|) for the SDM reflect the statistics of all possible 
consecutive two-click events separated by a time t. In 
order to judge about the correlations between two truly 
successive atoms, one should take the limit t — > in 
Eqs. H37|l and (|39|l . In this limit, two consecutive atoms 



interact with the cavity field and there is no time for 
decoherence to take place in between. Therefore, one ex- 
pects stronger atomic correlations as a consequence of 
the discussion done in Sec. Ill Bl 

In Fig. 01 we show the correlation function Geg as a 
function of the displacement We see that Geg exhibits 
stronger correlations, when compared with the conven- 
tional micromaser, despite the complete classicality of 
the SDM steady-state. This is a counter intuitive re- 
sult coming from the belief that stronger correlations 
should appear only when non-classical steady-states are 
involved. The stronger atomic correlations present in the 
SDM originate in the different nature of the unitary pro- 
cess. In the SDM, we rely on an interaction (see Sec. 
that naturally produces mesoscopic entangled atom-field 
states ( "Schrodinger cat states" ) , while the conventional 
micromaser uses the Jaynes-Cummings interaction, pro- 
ducing Rabi oscillations inside well defined atom-field 
subspaces, and essentially exchanging a single photon per 
cycle. 



IV. NUMERICAL RESULTS 

In order to describe the SDM dynamics we have at our 
disposal the analytical expression of the symmetrically 
ordered characteristic function x(/3,/3*,i), Ea. (PS|l . We 
can as well describe the time evolution of the Wigner 
function W{a,a* ,t), that is the Fourier transform of 
x{f3,f}*,t), which gives a picture in the phase space as- 
sociated to the cavity mode. We can further consider 
the time behavior of the density matrix {pF)m,n{t) — 
{m\pF{t)\n), derived from the master equation in Eq. 
when the field density operator is represented in the Fock 
basis. In this case, we can use quantum jumps techniques 
[l7j , already applied succesfuUy in several problems in the 
domain of CQED [Ulli. 

All these intertwined tools allow a consistent descrip- 
tion of SDM dynamics which complements the analytical 
results and unveils additional features. As an example, 
in Fig. [3 we show the transient behavior of the Wigner 
distribution M^(q;) of the SDM field. Starting from the 
Gaussian function of the vacuum state (Fig. |3Ji), after 
a time interval = 0.1 (Fig. ^p) we see the presence 
of two additional peaks, symmetrically placed along the 
imaginary axis. Here, we have chosen |^| — tt, a large 
enough value to see a peaked structure. At later times in 
the transition, Fig.|3fc,d), we see the onset of other peak 
pairs symmetrically placed on the imaginary axis, while 
the distribution lowers and broadens along that axis. The 
underlying physical mechanism for the generation of this 
dynamics was described in Sec. ^1 

At later times decoherence destroys such structures, 
and when the steady-state is reached (Fig. ^ the dis- 
tribution looks like the envelope of the multiply-peaked 
structure along the imaginary axis, while it preserves its 
initial minimum width along the real axis. FigureE|shows 
the density matrix of the steady-state SDM field. The di- 
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FIG. 5: Generation of multi-peaked distributions in the tran- 
sient cavity field dynamics; W{a,t) distributions at t = 
(top) and yt — 0.f,0.3, 0.5 (proceeding downward), obtained 
from the master equation I21II . 




FIG. 6: Steady-state cavity field after the transient evolution 
of Fig. El W{a) function at 'yt = 20. 
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FIG. 7: Density matrix of the steady-state cavity field. 



agonal elements provide the photon statistics, and their 
wide distribution confirms the predicted super-poissonian 
behavior. Furthermore, we note the presence of ofF- 
diagonal density matrix elements or coherences, which 
rule the phase and spectral properties of the field. 

Initially they are not excited in the cavity field, hence 
they are induced by the interaction with the strongly 
driven atoms. In order to investigate this effect, in Fig.|Sl 
we show the steady-state Pegg-Barnett ^ phase distri- 
bution 



i{m—7i)0 



(41) 



m,n— 



The SDM phase distribution shows a particular fea- 
ture, i.e., a narrow two-peaked structure centered on the 
values 9 = ±7r/2, which is already present in the tran- 
sient. This feature is explained by the SDM gain dynam- 
ics, Eq. H19|) . which shows that the resonant interaction 
of the cavity field with strongly driven atoms implies an 
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FIG. 8: Pegg-Barnett phase distribution of the steady-state 
cavity field. 

equal displacement of the field by the imaginary quanti- 
ties In turn, this effect is quite consistent with the 
narrow elongated form of the steady-state distribution 
along the imaginary axis in phase space, as well as with 
the vanishing of the steady-state field expectation value. 

V. CONCLUSIONS 

We have investigated the dynamics of a strongly-driven 
micromaser, based on the resonant interaction of one 
mode of a high-Q cavity with a poissonian low density 
beam of two-level atoms strongly driven by a resonant 
classical field. We have shown that this system pro- 
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